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We investigate the IceCube detection potential of very high energy neutrinos from blazars, for 
different classes of "hadronic" models, taking into account the limits imposed on the neutrino flux 
by the recent Fermi telescope observations. 

Assuming the observed 7-ray emission is produced by the decay of neutral pions from proton- 
proton interactions, the measurement of the time-averaged spectral characteristics of blazars in 
the GeV energy band imposes upper limits on the time-averaged neutrino flux. Comparing these 
upper limits to the 5cr discovery threshold of IceCube for different neutrino spectra and different 
source locations in the sky, we find that several BL Lacs with hard spectra in the GeV band are 
within the detection potential of IceCube. If the 7-ray emission is dominated by the neutral pion 
decay flux, none of the flat-spectrum radio quasars are detectable with IceCube. If the primary 
high energy proton spectrum is very hard and/or neutrinos are produced in proton-photon, rather 
than proton-proton reactions, the upper limit on the neutrino flux imposed by the measured 7-ray 
spectra is relaxed and gamma-ray observations impose only lower bounds on the neutrino flux. We 
investigate whether these lower bounds guarantee the detection of blazars with very hard neutrino 
spectra (spectral index F„ ~ 1), expected in the latter type model. 

We show that all the "hadronic" models of activity of blazars are falsifiable with IceCube. Fur- 
thermore, we show that models with 7-ray emission produced by the decay of neutral pions from 
proton-proton interactions can be readily distinguished from the models based on proton-gamma 
interactions and/or models predicting very hard high energy proton spectra via a study of the 
distribution of spectral indices of 7-ray spectra of sources detected with IceCube. 



I. INTRODUCTION 

The broadband, radio-to-7-ray spectra of blazars 
are most commonly interpreted within the frame- 
work of synchrotron-sclf-Compton (SSC)or synchrotron- 
external-Compton (SEC) models 0, H [| II I! ■ Within 
this framework, the observed 7-ray emission is gener- 
ated via inverse Compton (IC) scattering of low energy 
background photons. The background photons originate 
either from an external source (in SEC models) or are 
produced by the same high energy electrons via the syn- 
chrotron mechanism (SSC models). 

An alternative interpretation of the 7-ray emission is 
considered in models involving interactions of high energy 
protons H, 0, Hi- Several types of such "hadronic" inter- 
action models are considered in the literature. Different 
models adopt different assumptions about the dominant 
interaction channel of high energy protons and about the 
high energy proton spectrum. 

If high energy protons interacting with low energy pro- 
tons is the dominant channel (e.g. from the accretion 
flow onto a supermassive black hole), a significant con- 
tribution to the observed 7-ray flux could be secondary 
emission from the decay of neutral pions, 7r° — > 77. In- 
deed, if one assumes a "conventional" E~ 2 power-law 
type spectrum of high energy protons (see e.g. @) gen- 
erated for instance by proton shock acceleration in the 
blazar jet, one finds that the spectrum of 7-ray produced 
by the tt° decay should be distributed over a broad range 
of energies, starting from a fraction of the pp reaction en- 
ergy threshold, E pp ^hr = 2m w (l + m^/m p ) ~ 280 McV, 



where m p , are proton and pion masses, up to the 
high energy cutoff of the proton spectrum. In con- 
trast, if the high energy protons interact mostly with the 
background photons, the reaction threshold is kept high, 
£ P7ithr = {nipin^ + m 2 )/e ~ 6.8 x 10 16 [e/1 cVp 1 cV, 
where e is the energy of the background photon. Con- 
sequently, the observed 7-ray emission from blazars, de- 
tected at energies E 1 -C E Pit hr, cannot be produced by 
neutral pion decay. Instead, the observed 7-ray emission 
in this type of models is usually supposed to be produced 
by secondary electron-positron pairs generated in proton 
interactions and/or in electromagnetic cascade initiated 
by the proton interaction. 

A similar "proton-initiated cascade" origin of the ob- 
served 7-ray emission is assumed in models which con- 
sider the possibility of very hard spectra of high energy 
protons, or proton spectra with low energy cutoffs (see 
e.g. 0, EE 111 EE El EE El). Such hard proton spectra 
can be produced if protons are accelerated by large-scale 
electric fields induced close to the black hole horizon, 
rather than at shocks in the blazar jets. Arguments in 
favor of particle acceleration directly in the central en- 
gine of active galactic nuclei (AGN) are given by the 
recent observations of very fast variability of high en- 
ergy 7-ray emission from blazars and radio galaxies (see 

e.g.[IlE|E3)- 

A straightforward distinction between hadronic and 
leptonic models and between different types of hadronic 
models of emission from blazars would follow from the de- 
tection of high energy neutrinos, expected to accompany 
7-ray emission in hadronic models only. However, ex- 
perimental limitations have prevented positive outcomes 
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until recently. First, the sensitivity of existing high en- 
ergy neutrino detectors was not sufficient to detect neu- 
trino flux even from the brightest blazars. Search results 
for neutrino signal from astronomical point sources have 
until now only provided upper limits fl8l . Il9l . |20| usu- 
ally above model predictions. Next, the quality of the 
available 7-ray data was not sufficient to make firm pre- 
dictions of the time-averaged (over flaring and quiescent 
states, taking into account flaring "duty cycles") neu- 
trino fluxes from blazars expected in different hadronic 
emission models. 

With the newly operating IceCube [2l[ and Fermi [22j 
telescopes, a new sensitivity window on high energy neu- 
trinos is currently opening (see e.g. [23| where the Ice- 
Cube science potential has been generically described or 
[2I . regarding galactic sources) and the detection of 
quiescent (time-averaged) levels of 7-ray emission from 
a large number of blazars has become possible. In this 
paper we discuss the implications of the time-averaged 7- 
ray spectral measurements of a large number of blazars 
for the detection of neutrinos from this class of sources 
by IceCube. 

Under the assumption of an observed 7-ray emission in 
the GeV energy band dominated by neutral pion decay, 
the time averaged neutrino flux from blazars can be pre- 
dicted. Assuming a simple power-law type spectrum of 
the primary high energy protons, the measurement of the 
7T° decay flux and spectral index in the GeV band fully 
constrains the expected neutrino flux at higher energies. 
Using the first published Fermi catalogue of blazars and 
AGN (normalization and spectral characteristics) [2(| . 
the corresponding neutrino fluxes for each blazar are cal- 
culated. Once these neutrino fluxes convolved to the in- 
strumental response of IceCube for sources at different 
declinations, the expected number of neutrinos for each 
blazar is derived. This analysis leads us to conclude a de- 
tectable steady-state neutrino emission from a few (out 
of ~ 10 2 ) Fermi blazars with the hardest spectra in the 
GeV band. 

If the observed 7-ray emission is dominated by 
the synchrotron and/or inverse Compton and/or 
Bremsstrahlung emission from the secondary e + e~ pairs 
produced in proton-initiated particle cascades, the con- 
straints on the neutrino flux from the source are signifi- 
cantly weaker, since the measurement of the GeV band 
7-ray spectrum does not provide a direct constraint on 
the neutrino spectrum. We demonstrate this possibility 
by considering a model in which a nearly monocnergetic 
proton beam (composed of e.g. protons accelerated in 
a vacuum gap in a black hole magnetosphere) produces 
7-rays and neutrinos in proton-proton interactions (e.g. 
during the propagation through the accretion flow). 

By combining IceCube and Fermi data, we show that it 
should be possible to (a) confirm or disprove the hadronic 
models of activity of blazars and (b) firmly distinguish 
between models which predict very hard neutrino spec- 
tra (which could be either due to the low energy cut- 
offs in the proton spectra or due to the threshold effects 



of the proton-photon interactions) and models which as- 
sume relatively soft (i? _2 -like) neutrino spectra). This 
would proceed via a study of the distribution of neutrino- 
detected blazars in a 7-ray "Flux vs. photon index" di- 
agram in the GeV band: In the scenario with neutral 
pion decay 7-ray emission from proton-proton interac- 
tions all the blazars detected with IceCube should occupy 
the "high flux - hard spectrum" corner of the diagram, 
while they should be randomly distributed in the models 
based on pj interactions and/or proton spectra with low 
energy cutoffs. 

The plan of the paper is as follows: In Section |TT] 
we consider a specific hadronic model of blazar activ- 
ity which assumes 7-ray emission produced a power-law 
distribution of high energy protons (with a relatively soft 
spectra) in interaction with low energy protons from the 
accretion flow. We show that in this class of models the 
spectral characteristics of the GeV-band 7-ray emission 
are directly related to the spectral characteristics of neu- 
trino emission. In Section lTTTl we derive the time-averaged 
neutrino flux from blazars detected in the GeV energy 
band by the Fermi telescope and demonstrate that only 
blazars with the hardest GeV-band spectra can be ex- 
pected to be detectable as neutrino sources. 

Next, in Section IIVI we present a specific hadronic 
model which predicts a very hard neutrino spectrum. 
We demonstrate that in this class of models, measure- 
ments of the GeV band 7-ray spectra by Fermi impose 
only lower bounds on the neutrino flux, but do not con- 
strain the spectral characteristics of the neutrino signal. 
We discuss the numerous uncertainties affecting the neu- 
trino flux estimates. In Section[V]we investigate whether 
blazars with very hard neutrino spectra can be detected 
with IceCube. One immediate problem for the detection 
of sources with very hard spectra is related to the Earth 
opacity to neutrinos with energies above ~PeV, which 
dominate the source flux. The statistics of the neutrino 
signal at increasing source declinations is therefore re- 
duced. We then derive the IceCube spectral discovery 
curves for arbitrary source declinations. Finally, in Sec- 
tion [VI] we investigate whether the hadronic model(s) of 
blazar activity are falsifiable with IceCube and whether 
different types of hadronic models can be distinguished 
using a combination of IceCube and Fermi data. 



II. GEV 7-RAY EMISSION FROM NEUTRAL 
PION DECAY 

If the high energy protons are produced via shock ac- 
celeration, their spectrum is expected to be roughly a 
cutoff powerlaw with a spectral index r p ~ 2, dN p /dE ~ 
E~ Fp exp(—E/E PtCU t). If collision of high energy protons 
with low energy protons is the main interaction chan- 
nel (e.g. from the accretion flow onto supermassivc black 
hole, or with protons populating the blazar jet), both the 
7T° decay 7-ray and the neutrino spectra are cutoff power- 
laws with approximately the same indices, T u ~ T 7 ~ T p , 
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in a wide energy interval from ~ 0.1 GcV up to a cutoff 
at energies somewhat smaller than that of the primary 
proton spectrum E JiUtCUt ~ KE Pi cut with K ~ 0.1 [27| . 

Naively, the similarity between the 7-ray and neu- 
trino spectra should enable to predict the blazar neutrino 
emission characteristics if the 7-ray emission spectrum is 
known. The IceCube telescope being most sensitive to 
neutrinos with energies above ~ 1 TeV, the knowledge 
of the 7-ray source spectrum in this energy band is most 
important to estimate the detectability of the neutrino 
signal. However, in the case of extragalactic sources, the 
7-ray flux at energies above ~ 10 GeV - 1 TeV is atten- 
uated by pair production on the background radiation 
fields (collectively known as Extragalactic Background 
Light, EBL). This degrades the information carried by 
the shape and overall flux of the ir° decay spectrum of 
the source. This, in turn, prevents a direct assessment 
of source detectability with a neutrino telescope. Even 
for relatively low redshift sources, such as TeV 7-ray - 
loud blazars, measurements of source spectra above 100 
GeV do not provide information about the source neu- 
trino spectra. This is illustrated in Fig. [TJ in which we 
show a model fit for the GeV- TeV band emission spec- 
trum of the blazar PG 1553+11, recently measured by 
the Fermi, HESS and MAGIC telescopes. The TeV band 
observations result only in a lower bound on (rather than 
measurement of) the source 7-ray flux (see e.g. (28j). 

Contrary to TeV 7-rays, 7-rays in the GeV energy band 
do not suffer from the absorption on the EBL. Measure- 
ments of the spectral characteristics in this energy band 
provide, therefore, information about the intrinsic spec- 
trum of the blazar. Since the ir° decay spectrum is ex- 
pected to extend down to sub- GeV energies in models 
based on proton-proton interaction, a measurement of 
the spectral characteristics in the GeV energy band could 
be used to constrain the ir° decay spectrum. Assuming 
a power-law primary proton spectrum, the intrinsic 7- 
ray spectrum of the source in the TeV energy range can 
be reconstructed and the neutrino flux from the measure- 
ments of the flux and spectral index of the 7-ray emission 
in the GeV band can be estimated. 

This possibility is illustrated in Fig. [1] in which we 
show the results of a numerical calculation of the 7-ray 
and neutrino spectra produced in proton-proton inter- 
action of high energy protons with a power-law spec- 
trum with index T p = 1.8 and a high energy cutoff at 
E p> cut = 10 18 eV, with low energy protons from the (ra- 
diativcly inefficient) accretion flow onto a 10 7 M Q black 
hole. Our model is identical to the model of high en- 
ergy emission from the cascade developing at the base 
of the jet of a TeV blazar, PKS 2155-304, developed in 
the Ref. Protons are initially accelerated close to 

the black hole horizon (at the distance D ~ 10 12 cm) 
and escape to large distances (our calculation stops at 
D = 10 16 cm). Interaction with protons of the accretion 
flow (see [Iff for the detailed description of the geome- 
try of the radiativcly inefficient accretion flow) initiate 
a particle cascade which accompanies the escape of the 
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FIG. 1: Spectra of 7-ray and neutrino emission of blazar 
PG 1553+11 expected in a model in which the GeV emis- 
sion is dominated by the 7-rays from neutral pion decay. The 
thin solid, short dashed and dotted lines show the contribu- 
tions from synchrotron, inverse Compton and Bremsstrahlung 
emission from the secondary e + e~ pairs. The two curves 
marked IceCube show the discovery potential of IceCube to 
the power-law like neutrino spectra with different spectral in- 
dices, for different source declinations: 15° (thick solid red 
curve) and 75° (thick dashed red curve), derived in Section 



high energy protons. The cascade leads to the gener- 
ation of neutrino and 7-ray emissions, calculated using 
the analytical approximations for the production spectra 
of secondary particles in pp interaction derived by Kclner 
et al. 0. 

From Fig. [1] one can see that the GeV band emission 
is dominated by (or, at least, has a significant contri- 
bution from) the 7r decay emission component, which 
has a spectral energy distribution very similar to the one 
from neutrinos. This enables an estimate of the charac- 
teristics of neutrino spectrum from the GeV-band 7-ray 
spectrum. 

One should note that, apart from the tt° decay emis- 
sion, the 7-ray spectrum in the GeV band could have 
contributions from synchrotron, inverse Compton and 
Bremsstrahlung emission produced by the secondary 
e + e~ pairs. Such e + e~ pairs are produced by charged 
pion decay and/or in result of the development of elec- 
tromagnetic cascade inside the blazar jet. In the pure 
hadronic model of blazar activity, where all the high en- 
ergy electrons and positrons are supposed to be produced 
in this way, the power injected in electrons/positrons 
from the tt^ decay is at most comparable to the power 
injected into ir° decay 7-rays. If the primary proton spec- 
trum is relatively soft, r p ~ 2, the power transferred to 
the e + e~ pairs by the electromagnetic cascade is also 
at most comparable to the initial ir° decay 7-ray power. 
This means that in models based on proton-proton inter- 
action with soft high energy proton spectra, a significant 
contribution in the GeV band 7-ray emission is given by 
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FIG. 2: Distribution of Fermi AGN in "flux vs. spectral in- 
dex" parameter space (from Abdo et al. HH). Grey curves 
show the 5<r "discovery" threshold for detection of neutrino 
signal from Fermi blazars by IceCube in the model in which 7- 
ray emission is dominated by the signal from neutral pion de- 
cay. Different curves correspond to different assumed source 
declinations, from 8 — 5° (upper curve) to 5 = 85° (lower 
curve) with a AS — 5° step. Below the 3-yr IceCube discov- 
ery curves are indicated the 1-yr ANTARES 90% C.L upper 
limit curves, extracted and extrapolated from [33] (see text). 
Blue filled circles correspond to Fermi blazars considered as 
"neutrino source candidates" in the model with neutral pion 
decay 7-ray emission. 

the 7T decay component, in spite of the possible presence 
of synchrotron, inverse Compton and/or Brcmstrahlung 
emission components. This ensures the possibility to es- 
timate the neutrino signal from the detected GeV 7-ray 
signal. 

III. DETECTABILITY OF QUIESCENT 
NEUTRINO EMISSION FROM BLAZARS IN 
THE LIGHT OF FERMI 

Measurements of the spectra of quiescent emission 
from blazars at energies below 10 GeV became possible 
recently with the start of operation of the Fermi tele- 
scope. Measurements of the spectral parameters of the 
7-ray emission from blazar, reported in po| enables, for 
the first time, to make predictions of the time-averaged 
neutrino flux within the hadronic model of activity, un- 
der the assumption that the GeV band flux is dominated 
by the ir° decay contribution. 

The resulting estimates of the number of neutrinos 
from different blazars, given three years of IceCube ex- 
posure (in the 80 string baseline configuration) [2l| are 
shown in Table U (only blazars with an expected number 
of neutrinos close or above the 5a " discovery" threshold 
are listed). 

The 5cr "discovery" threshold is calculated using the 
declination-dependent atmospheric neutrino background 
parametrization from [29j (which includes the "prompt" 



Pfi flux). The "prompt" component is as yet uncer- 
tain and may significantly contribute to the overall at- 
mospheric neutrino spectrum at energies above E ~ 
10 5 ~ 6 GeV. 

The diffuse neutrino emission from the interaction of 
galactic cosmic rays with the interstellar medium, es- 
timated for instance in Ref. [3(|, has been neglected: 
While expected to be relatively important for sources 
at low galactic latitudes 6^0°, but not exceeding the 
"prompt" atmospheric neutrino component, it is largely 
sub-dominant at large galactic latitudes. 

The diffuse neutrino emission from unresolved blazars 
was neglected as well. Within the proton-proton in- 
teraction models considered above and based on early 
measurements of the EGRET extragalactic diffuse 7-ray 
background by (3lj . the diffuse extragalactic neutrino 
flux was estimated to eventually dominate the atmo- 
spheric neutrino flux above E ~ 10 5 GeV (32j. How- 
ever, revised measurements of the extragalactic diffuse 
7-ray background with EGRET indicate a lower normal- 
ization and a softer background spectrum (33j . thus ac- 
cordingly reducing this neutrino background component 
from a rcscaling of the estimates in Ref. [32| to a level 
well below the atmospheric neutrino flux at ~ 10 5 GeV. 

The uncertainty of the diffuse neutrino background 
level at energies above E ~ 10 5 GeV may slightly mod- 
ify the results of Table [I] especially the numbers of 
neutrino-induced muon events with detector threshold 
energy E™ in - 2 = 10 4 5 GeV. On the other hand, the 
numbers of neutrino-induced muon events with detector 
threshold energy E™ 1 ™* 1 = 10 3 - 5 GeV remain almost un- 
affected, the dominant diffuse background source being 
the conventional atmospheric neutrino component. 

To estimate the parameters of the neurino spectrum 
(normalization <£>„, and spectral index T„) from the ob- 
served parameters of the 7-ray spectrum (<j> 7 ,r 7 ), we 
produce a "lookup table" (r v , $„) vs. (r 7 ,$ 7 ) using 
the neutrino and 7-ray production spectra for different 
powerlaw distributions of the prima ry p rotons. We use 
the approximations of Kelner et al. [34| for the calcula- 
tion of the production spectra and assume a fully oscil- 
lated neutrino flux (i.e. the normalization of the muon 
neutrino spectrum is 1/3 of the normalization of the 
primary total neutrino spectrum). The number of de- 
tected neutrinos from the source (and the contamina- 
tion by the atmospheric neutrino background) depends 
on the low muon energy threshold chosen for the analy- 
sis. In Table [H we show the number of neutrinos from 
the source for two different choices of the reconstructed 
muon energy thresholds, E™ 171 ' 1 = 10 3 5 GeV ~ 3TeV 
and E™ n - 2 = 10 4 5 GeV ~ 30TeV. In analyzing the 
data, the best choice for the reconstructed muon energy 
threshold depends on the neutrino spectral index, V v . 
This issue is discussed in greater details in Section IVl 

Out of 104 bright blazars detected by Fermi during 
the first three months of operation [26[, only a few are 
expected to produce sufficiently high neutrino flux for 
a detection within one year of IceCube operation. All 
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the neutrino candidate blazars are BL Lac type objects. 
None of the bright and powerful Flat Spectrum Radio 
Quasars (FSRQ), such as 3C 279, 3C 273 or 3C 454.3, are 
expected to produce a neutrino flux above the IeeCube 
detection threshold. 

This is also clear from Fig. [2]where we show the thresh- 
olds of detectability of neutrino emission from Fermi 
blazars, calculated by normalizing the neutrino flux from 
the Fermi 7-ray spectra. In principle, the expected num- 
ber of high energy neutrinos from the source depends not 
only on the intrinsic source spectrum, but also on the po- 
sition of the source in the sky (given the specific location 
of IeeCube at the South Pole, on the source declination). 
Not only the number of source neutrinos, but also the 
atmospheric neutrino flux, which constitutes the back- 
ground for the source signal, depend on the declination. 
Thus, it follows a dependence of the 5a source detection 
threshold on the declination. Fig. [2] shows that the best 
candidates for a detection with IeeCube (within proton- 
proton models with soft high energy proton spectra) are 
the low declination BL Lacs with hard GeV band 7-ray 
spectra. All the bright FSRQ detected by Fermi are ex- 
pected to produce neutrino counts below the detection 
threshold of IeeCube, independent of their declination. 

To understand why the bright FSRQ are not detectable 
by IeeCube, one has to take into account that the pre- 
dictions for the neutrino flux depend not only on the 
normalization of the 7-ray spectrum at GeV energies, 
but on its hardness as well. IeeCube is most sensitive 
to neutrinos with energies above ~ 1 TeV. The harder 
the spectrum of the source in the GeV energy band, the 
higher the expected flux in the TeV-PeV band. It turns 
out that in the GeV band, FSRQ have spectra on average 
softer than BL Lacs and thus explains why the less lumi- 
nous (in the GeV band) BL Lacs are expected to produce 
higher neutrino fluxes in the TeV-PeV band than FSRQs. 
These estimates rely however on an initially foreseen 80 
string IeeCube configuration and enhanced sensitivity to 
lower energy in the new 86 string baseline configuration 
has not been considered yet [H[ . 

It is worth mentioning that most of the best candidate 
sources listed in Table U are TeV 7-ray loud BL Lacs. 10 
out of 16 sources have already been reported to be TeV 
emitters. The remaining 6 have either large or unknown 
redshift which explains the difficulty of their detection 
above 100 GeV. It is unlikely that the TeV blazar neu- 
trino emission be produced in pry interaction, because a 
high efficiency of the pj interaction (which has a cross 
section a pl ~ 10~ 28 cm 2 ) would automatically imply a 
strong absorption of the TeV 7-ray s via pair production 
on the soft photon background in the source (the cross 
section for this process is a 11 ~ 10~ 25 cm 2 , i.e. three 
orders of magnitude larger than the pry cross section) 
[l3l [36|. The only possibility for neutrino production 
in TeV blazars is via proton-proton interactions. 

The detection potential of ANTARES for the strongest 
neutrino emitters can also be roughly estimated using the 
one year E~ 2 90% C.L upper limit curves found in Ref. 



[13 (E 2 d<5> u J&E u « 4 - 10 x 10~ 8 GeV cm" 2 s" 1 over 
the declination range —90° < 6 < 30°), assuming a be- 
havior similar to IeeCube for hardening spectra (i.e. by 
extrapolating parallel to IeeCube 5a discovery curves). 
The expected signal escapes ANTARES discovery poten- 
tial for all sources but for PG 1553+11, which could be 
in the reach of ANTARES after one year. Note how- 
ever that this assumption may not hold, see Appendix 
lAl given the distinct detector response of ANTARES and 
IeeCube. 



IV. HADRONIC MODELS WITH VERY HARD 
NEUTRINO SPECTRA 

The above IeeCube neutrino count estimates from 
blazars show that within the hadronic models predicting 
relatively soft neutrino spectra (E~ 2 like), < 10 7-ray 
loud blazars with hardest GeV band 7-ray spectra can 
be detected within a several year- long neutrino data set. 

Although experimental results in the literature usu- 
ally consider benchmark power-law type neutrino spec- 
tra dN v /dE (x E~ T " with r„ ~ 2 most existing hadronic 
models of blazar activity predict neutrino spectra which 
significantly differ from a E~ 2 powerlaw. This is often 
related to the existence of a low energy cutoff in the pro- 
ton spectra. As it is mentioned in the Introduction, these 
low energy cutoffs can be either intrinsic cutoffs in the 
proton spectra (e.g. in the models with proton accelera- 
tion in the blazar central engine) or can arise due to the 
existence of a threshold energy of proton-photon interac- 
tions. 

Since the low energy cutoff or the photo-pion thresh- 
old energy, considered in most of the hadronic mod- 
els, is much higher than ~ 1 PeV, the contribution of 
the 7T° decay component to the 7-ray emission spec- 
trum in the GeV - TeV energy band is small, com- 
pared to the synchrotron and/or inverse Compton and/or 
Brcmsstrahlung emission from the secondary e + e~ pairs 
produced in the decay of the charged pions and/or in elec- 
tromagnetic cascade. Therefore, the estimates of the ex- 
pected spectral characteristics of neutrino emission based 
on the spectral characteristics of the GeV band 7-ray 
emission, presented above, do not work. 

To demonstrate this, we show in Fig. [3] the result of 
calculation of the 7-ray and neutrino spectra in the same 
model as in Fig. [TJ but assuming a very hard high en- 
ergy proton injection spectrum, with T p ~ (i.e. most 
of the power of the primary proton beam is carried by 
protons with energies close to the assumed cutoff energy 
E Pt cut = 10 18 eV). One can clearly see that a satisfactory 
fit of the GeV- TeV 7-ray spectrum is achieved with the 
dominant contribution from synchrotron emission from 
the secondary electrons. No information on the tt° decay 
spectrum, and hence no information about the neutrino 
spectrum can be extracted from the available 7-ray data. 

The hypothesis of hard neutrino spectra opens a pos- 
sibility of strong neutrino emission, even in the quiescent 
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Name 


Fermi Name z 


r 7 


F 




(n s , Tib, U5a-) E rrLin,2 




Blazars above 5a 


detection threshold 




B3 0133+388 


OFGL J0136.6+3903 (?) 


1.65 


1.8 


(14, 0.36, 7) 


(4.6, 0.0079, 3) 


1ES 1011+496 OFGL J1015. 2+4927 0.212 


1.73 


4.9 


(12, 0.31, 7) 


(3.1, 0.0067, 3) 


Mrk 421 


OFGL J1104.5+3811 0.0300 


1.77 


15.3 


(27, 0.36, 7) 


(7.7, 0.0080, 3) 


PKS 1424+240 


OFGL J1427.1+2347 (?) 


1.80 


6.2 


(11, 0.53, 8) 


(3.6, 0.013, 3) 


PG 1553+11 


OFGL J1555.8+1110 0.36(?) 


1.70 


8.0 


(71, 0.95, 10) 


(31, 0.024, 4) 


Mrk 501 


OFGL J1653.9+3946 0.0337 


1.70 


3.1 


(13, 0.35, 7) 


(3.9, 0.0077, 3) 


PKS 1717+177 


OFGL J1719.3+1746 (?) 


1.84 


6.9 


(7.9, 0.67, 8) 


(2.8, 0.017, 4) 




Blazars close to 5a detection 


threshold 




3C 66A 


OFGL J0222.6+4302 0.444 


1.97 


25.9 


(3.3, 0.34, 7) 


(0.74, 0.0076, 3) 


AO 235+164 


OFGL J0238.6+1636 0.94 


2.05 


72.6 


(5.2, 0.70, 9) 


(1.5, 0.018, 4) 


1ES 0502+675 


OFGL J0507.9+6739 0.416 


1.67 


1.7 


(4.8, 0.24, 6) 


(0.91, 0.0038, 3) 


B 1218+30 


OFGL J1218. 0+3006 0.182 


1.89 


9.7 


(4.3, 0.43, 7) 


(1.2, 0.010, 3) 


W Com 


OFGL J1221. 7+2814 0.102 


1.93 


8.3 


(2.3, 0.46, 7) 


(0.63, 0.011, 3) 


1ES 1959+650 OFGL J2000. 2+6505 0.047 


1.86 


4.2 


(1.3, 0.24, 6) 


(0.20, 0.0041, 3) 




Neutrino candidate blazars in 


Southern hemisphere 




KUV 0311-1938 


OFGL J0033.6-1921 0.610 


1.70 


1.6 






PKS 2005-489 


OFGL J2009.4-4850 0.071 


1.85 


2.9 






PKS 2155-304 


OFGL J2158.8-3014 0.116 


1.85 


18.1 







TABLE I: List of blazars with potentially detectable quiescent neutrino emission, in the model with 7-ray flux in the Fermi 
band dominated by neutral pion decay. Bold names: blazars detected in the TeV energy band. Blue (Grey) names: blazars 
in the Southern hemisphere, not accessible to IceCube. The two last columns list the expected number of detected signal and 
atmospheric neutrinos, n a and rib, together with the estimated 5a discovery threshold, n^, a for two choices of the reconstructed 
low muon energy threshold. E™ inA and E™ in < 2 (see text). 



T 



-8 - 




6 7 8 9 10 11 12 13 1 A 15 16 17 18 19 

log E [eV] 

FIG. 3: Spectra of 7-ray and neutrino emission of blazar PG 
1553+11 expected in a model with GeV emission dominated 
by 7-rays from secondary e + e~ pairs. The notations are the 
same as in Fig. [T] 



state, from a large number of blazars, including both BL 
Lacs and FSRQs. 

It was argued in Ref. [38| that even if the GeV - 
TeV band 7-ray emission is dominated by the flux from 
electron-positron pairs from tt^ decay, constraints on the 
neutrino luminosity of the source can still be derived from 



the spectra of the broad band 7-ray emission. Namely, 
a measurement of the overall bolomctric luminosity of a 
Galactic 7-ray source imposes a restriction on its neu- 
trino luminosity (but not on the spectral index or cutoff 
energy of the neutrino spectrum). 

In the case of extragalactic sources, measurements of 
the electromagnetic luminosity of the source also imposes 
a restriction on its neutrino luminosity. However, the re- 
lation between the observed 7-ray power and morphology 
of an extragalactic source and its neutrino luminosity is 
more complicated than in the case of a Galactic source. 

The problem being that it is possible that even for 
sources at small rcdshifts, most of the 7-ray power ini- 
tially emitted in the energy band above 0.1-1 TeV is 
absorbed during propagation through the intergalactic 
medium, due to pair production on the EBL. This means 
that the primary source power produced in result of pro- 
ton interaction is not known. At the same time, to ob- 
tain a correct estimate of the neutrino luminosity, one 
has to be able to estimate the primary source power at 
energies above 0.1-1 TeV. This could, in principle, be 
done, because the absorbed power is initially converted 
into e + e _ pairs and then back to the 7-ray s, because the 
e + e~ pairs loose energy via inverse Compton scattering 
on the CMB photons. Measuring the flux of the inverse 
Compton emission from the 7-ray-induccd cascade e + e~ 
pairs one would be able to measure the flux of the initial 
7-ray s and in this way to measure the 7-ray luminosity 
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of the primary source. 

The method of measurement of the flux from the 7- 
ray-induced cascade e + e~ pairs heavily depends on the 
strength of the extragalactic magnetic fields. If the ex- 
tragalactic magnetic fields are strong. B 3> 10~ 12 G, the 
trajectories of the e + e~ pairs produced in the process of 
absorption of the 7-rays on the EBL would strongly de- 
viate before they transfer energy through inverse Comp- 
ton scattering to the CMB photons [39(. The emission 
from the 7-ray-induced cascade should therefore be de- 
tectable as an extended 7-ray halo around the primary 
point source. The size of the halo is determined by the 
mean free path of the 7-ray s through the EBL and thus 
depends on the primary 7-ray energies. The luminosity 
of the halo is equal to the 7-ray luminosity of the pri- 
mary source in the TeV - EeV energy band, averaged 
over 10 3 — 10 6 yr time scale. 

If the extragalactic magnetic fields are somewhat 
weaker, 10" 16 < B < 1(T 12 G, the deflection of e+e" 
pairs over the radiative cooling time would be moder- 
ate, but still the signal from the 7-ray induced cascade 
in the intergalactic medium could also be detected as an 
extended emission region around an extragalactic point 
source [4(| EH , with the source size having a somewhat 
different dependence on the energy than for halos dis- 
cussed by Aharonian et al. . 

Finally, if the magnetic fields in the intergalactic 
medium are B <C 10~ 16 G, the extended emission from 
the 7-ray-induced cascade in the intergalactic space is in- 
distinguishable from the point source emission, because 
the size of the extended source becomes smaller than 
the point spread function of the current generation 7- 
ray telescopes [42[. I n this case a measurement of the 
point source flux provides a good estimate of the total 
electromagnetic and neutrino power of the source. 

One should note, however, that at present the strength 
of the extragalactic magnetic fields is largely unknown 
and the expected extended 7-ray emission from the ex- 
tragalactic sources is not detected. This means that the 
luminosity of the GcV - TeV 7-ray loud blazars at ener- 
gies above several TeV is largely unconstrained and, as 
a consequence, the neutrino luminosity of the blazars is 
not constrained by the 7-ray data. 

The overall luminosity of the source, generated by the 
interaction of high energy protons is larger or equal to 
the power contained in electromagnetic cascades, which 
transfer the initial source power to the sub- TeV energy 
band. Since the total electromagnetic luminosity of the 
source is comparable to its neutrino luminosity, measure- 
ment of the source flux in the sub- TeV energy band pro- 
vides a lower bound on the neutrino luminosity of the 
source. Taking this fact into account we show the cal- 
culated neutrino flux in Fig. [3l as a lower limit on the 
expected neutrino flux, obtained under the assumption of 
an observed 7-ray luminosity of the source comparable to 
its intrinsic 7-ray luminosity. 

From the discussion of this section it is clear that the 
search for very hard neutrino sources is extremely im- 



portant in the context of the study of blazar activity 
mechanisms. It is not clear a priori, if neutrino sources 
with very hard spectra, with most of the energy output in 
the > 10 PeV energy band, are accessible for study with 
IceCube. The problem being that at energies > 10 PeV 
the Earth is not transparent to neutrinos. Thus most of 
the signal from the source is absorbed in the Earth and 
does not reach the detector. To the best of our knowl- 
edge, no comprehensive study of sensitivity of IceCube to 
sources with very hard neutrino spectra was published up 
to now. Taking this into account, we investigate in the 
following section the potential of IceCube for the detec- 
tion of neutrino sources with very hard spectra extending 
to the multi-PeV - EeV energies. 



V. ICECUBE SENSITIVITY FOR 
ASTRONOMICAL SOURCES WITH HARD 
NEUTRINO SPECTRA 

The sensitivity of a neutrino telescope to a flux model 
can be derived from the knowledge of the neutrino ef- 
fective area of the telescope at the source declination 
(from the available information in publications, usually 
a declination-averaged effective area) and the angular & 
energy resolution of the instrument together with the as- 
sumption of an atmospheric neutrino flux. Nevertheless, 
this compact approach suffers limitations owing to the 
fact that the knowledge of the neutrino effective area 
alone does not permit the derivation of the optimal sen- 
sitivity to an arbitrary neutrino flux shape, as detailed 
information on the detector response is hidden behind 
the neutrino effective area concept. 

A refined method to assess the potential of a neutrino 
telescope to an arbitrary flux and arbitrary declination 
is based on the detector muon effective area and differ- 
ential muon flux at the detector, which can both be re- 
spectively extracted from the averaged neutrino effective 
area and from the neutrino flux spectrum propagated 
through the Earth following the methodology introduced 
by Neronov & Ribordy [381 ] . From the interplay of neu- 
trino signal contamination by atmospheric neutrinos and 
neutrino absorption, the optimal reconstructed muon en- 
ergy threshold for highest sensitivity for arbitrary spec- 
tral index and declination can be calculated: while for 
hard spectra it is advantageous to suppress low energy 
events dominated by atmospheric neutrinos in order to 
increase the discovery potential, with larger declinations 
the increasing neutrino absorption may play an impor- 
tant role above 10 - 100 TeV on defining the optimal 
reconstructed muon energy threshold. 

Based on this method, not only the best sensitivity 
can be derived but the possibility of detailed study of 
the source spectrum is also opened, provided sufficient 
statistics avoiding the unfolding of the measured muon 
spectrum. We apply the recipee from Neronov & Ri- 
bordy [38] for the calculation of the event rate in terms 
of reconstructed muon instead of neutrino energies. This 
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approach presents several advantages including the possi- 
bility of exploiting some of the physical observables (such 
as reconstructed energy) for the optimization of IceCube 
performance for the detection of sources with different 
spectral characteristics. 



A. IceCube "Spectral Discovery Curves" 

A major challenge for the detection of astronomical 
neutrino sources with IceCube is the separation of the 
source signal from the atmospheric neutrino background. 
If the spectrum of the source is harder than that of the 
atmospheric neutrino background, much of this back- 
ground could be suppressed by selecting only higher- 
energy events, i.e. adjusting the reconstructed low muon 
energy threshold in the data analysis event selection. The 
optimal reconstructed low muon energy threshold, which 
maximizes the signal-to-noise ratio, depends, in general, 
on the shape (e.g. on the spectral index, for a power 
law neutrino energy distribution) of the spectrum of the 
source neutrinos. The dependence of the optimal anal- 
ysis parameters on the spectral properties of the source 
leads to a dependence of the detector sensitivity on the 
assumptions about the spectral properties of the source. 
As a result, the knowledge of the IceCube discovery limit 
(= normalization of the neutrino spectrum at a given 
reference energy) for a "reference" source spectral index 
does not allow to estimate the sensitivity limit for a dif- 
ferent source spectrum. 

In Appendix [A] wc demonstrate a convenient way of 
presenting the IceCube discovery limits for arbitrary as- 
sumptions about the neutrino spectrum slope of the 
source. The "spectral discovery curves" shown Fig. @]for 
various source declination bear the following meaning: A 
neutrino flux with spectral index T„ (corresponding to a 
straight line in the \og(E) vs. E 2 log( Flux) representa- 
tion) tangent to one of the spectral discovery curve is at 
the limit of discovery. This representation of the IceCube 
discovery limit(s) is useful, because it enables 

• The estimate of the IceCube sensitivity for arbi- 
trary assumptions about the slope of the neutrino 
spectrum (i.e. the normalization of the minimal 
detectable spectrum at a given reference energy) 

• The estimate of the neutrino energy range con- 
tributing most significantly to the source signal 
(the energy at which the minimal detectable source 
spectrum touches the discovery curve). 

The dependence of the discovery curves, Fig. [4l on 
the source location in the sky is reduced to the source 
declination 8, given the privileged location of IceCube 
at the South Pole. We notice a weak dependence of the 
discovery limit on the source declination in the case of 
soft neutrino sources (r„ ^> 2): Most events contribut- 
ing to the source signal have relatively low energies, at 
which Earth is transparent for neutrinos. For sources 




log(E /GeV) 



FIG. 4: IceCube point source reach after three years of ex- 
posure. The tangent to one of these curves intersecting the 
y-axis at 1 GeV defines the muon neutrino flux normalization 
for the corresponding declination. 

with hard spectra, <C 2, the IceCube discovery po- 
tential is strongly affected by the increase of the source 
declination. The rise of the Earth absorption probabil- 
ity with energy obscures the source signal, contributed 
mostly by high energy neutrinos. The effect is most dra- 
matic for the hardest neutrino spectra, ~ 1. In this 
case, the minimal detectable flux for the sources at the 
declination 8 = 75° is about 1.5 orders of magnitude 
higher than from sources at 8 = 15°. 



VI. DISCUSSION: ICECUBE POTENTIAL FOR 
THE STUDY OF THE ORIGIN OF HIGH 
ENERGY BLAZAR ACTIVITY 

The IceCube sensitivity curves, shown in Fig. [¥] can 
be directly applied to estimate the detectability of neu- 
trino emission from Fermi blazars both in the model 
with 7T° decay origin of the 7-ray spectrum and with 
the 7-ray emission from the secondary e + e~ pairs. We 
demonstrate this on the example of the model calcula- 
tions for the TcV 7-ray loud blazar PC 1553+11. The 
red thick curves in Figs. [I] and [3j marked "IceCube", 
show the IceCube discovery curves for the source decli- 
nations 8 = 15° (close to the actual declination of PC 
1553+11) and 8 = 75°. The expected neutrino spectrum 
(thick green long-dashed line) clearly intersects the Ice- 
Cube discovery curves in both figures. Therefore, if the 
source activity is determined solely by acceleration and 
subsequent interactions of protons, the source should be 
detectable independent of the assumptions made about 
the origin of the GeV band 7-ray flux: The detection or 
non-detection of several blazars with IceCube will con- 
firm or rule out the so-called "hadronic" models of blazar 
activity. 

Note that, although the expected neutrino signal from 
PG 1553+11 in Figs. [1] and [3] were calculated within 
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the framework of a particular model of high energy pro- 
ton interaction in the source (interaction with low energy 
protons from a radiatively inefficient accretion flow) de- 
veloped in Ref . [l5[ , the numerical estimates of the over- 
all flux and spectral characteristics of the neutrino signal 
from the sources are, in a sense, model independent. 

Indeed, on the one hand, in all hadronic models of 
blazar activity, which assume the high energy 7-ray emis- 
sion dominated by the emission from the neutral pion 
decay, the normalization and the spectral characteristics 
of the neutrino flux are constrained by the measurements 
of the 7-ray flux in the GeV band in a similar way. 

On the other hand, in all hadronic models which 
assume that the GeV band 7-ray emission is domi- 
nated by synchrotron and/or inverse Compton and/or 
Bremsstrahlung emission from the e + e~ pairs produced 
in the charged pion decay, the measurement of the GeV 
band 7-ray flux imposes a lower bound on the neutrino 
luminosity of the source (the neutrino luminosity should 
be larger or equal to the measured 7-ray luminosity) . The 
shape of the neutrino spectrum at energies much below 
the threshold or low energy proton spectrum cutoff does 
not depend on the details of the proton spectrum itself: 
It is determined by the energy dependence of the dif- 
ferential cross section of proton interaction and charged 
pion decay spectrum. This means that in any model 
with 7-ray emission from secondary e + e~ pairs the ex- 
pected neutrino flux is at least as high as the one shown 
in Fig. [3] and therefore detectable with IceCube, pro- 
vided a low energy cutoff of the proton spectrum or the 
threshold energy of proton-photon interaction not higher 
than - 10 18 cV (see Fig. |3J). 

To summarize, we have explored the IceCube potential 
introducing a generic method to derive power-law type 
spectra discovery thresholds for an arbitrary neutrino 
spectral index (1 < r„ < 4) and for an arbitrary dec- 
lination. In the light of Fermi, the IceCube potential in 
the search for blazar neutrino point sources is promising 
even for quite different types of hadronic models of activ- 
ity: After three years of exposure, several sources should 
be revealed independently of the details of the hadronic 
models. Detection or non detection of the TeV sources 
PG 1553+11 and Mrk421 with the data collected by the 
IceCube collaboration between 2007 and 2009 would al- 
ready set strong constraints on the model with the soft 
spectrum of proton injection. 

Moreover, we established that IceCube will provide a 
definite answer regarding the validity of the hadronic 
class of models of activity, ruling them out in the case 
of no discovery. 

In contrast, the first high energy neutrino source detec- 
tions with IceCube would not only open a new observa- 
tional window, IceCube would delineate the two classes 
of hadronic models of activity as well, from a concurrent 
study of the Fermi photon and IceCube neutrino indices: 
While detected blazars conforming to models predicting 
soft neutrino spectra will lie in the lower part of the dia- 
gram Fig. [5] (delimited by the IceCube curves) , the ones 



conforming to models with hard neutrino spectra will be 
distributed over a large range of photon indices T 7 . The 
detection of blazars with IceCube would therefore truly 
marks the birth of multimessenger high energy neutrino 
astronomy. 



Acknowledgment 

The authors thanks L. Demirors and the anonymous 
referee for valuable comments. M. Ribordy is sup- 
ported by the Swiss National Research Foundation (grant 
PP002-114800). 



APPENDIX A: DERIVATION OF THE ICECUBE 
SPECTRAL DISCOVERY CURVES 

The neutrino point source strength normalization for 
discovery as a function of the neutrino source spectral in- 
dex can be derived from the calculations of the number of 
atmospheric and extra-terrestrial neutrino-induced muon 
events. These limits turn out to critically depend on the 
angular and energy resolution functions of the apparatus 
due to the diffuse nature of the soft atmospheric neutrino 
background. The relative importance of the resolution 
functions depend strongly on the spectral shape of the 
point source neutrino flux. The sensitivity to hard spec- 
tra with cutoff energy above ~ PeV will predominantly 
depend on the energy resolution while soft spectra will be 
probed best with instruments showing improved angular 
resolutions (Existing instruments deployed either in ice 
or water do not combine optimal resolution functions, 
which largely depend on the medium properties. This 
holds to the Cherenkov photon "thermalization" power 
of the media, determined by the relationship between 
Cherenkov light absorption and scattering lengths: In 
ice, A abs > A scat , in water A abs — A scat . Therefore, while 
instruments using ice as the Cherenkov photon transport 
medium show improved calorimetric properties, instru- 
ments deployed in water have better pointing capabilities 
for opposite reasons) . 

Our simplified calculation for the required number 
of observed events pointing back to a source for a 5a 
C.L. discovery is based on Poisson probability, requiring 
Cp(n b s — 1, Tib) > erf(5/v / 2) where n s (rib) is the number 
of signal (atmospheric) neutrino events, n Q b s = n s + 77.5 
and Cp is the Poisson cumulative distribution. 

For neutrino fluxes with spectral indices between 1 < 

r„ < 4, 



dE v 



E„ 



$ x (-^_) 

K GeV> 



(Al) 



we calculate the "5c discovery" flux normalization for 
a series of reconstructed muon energy thresholds, E™ m . 
Fig. [5] illustrates the procedure leading to the derivation 
of the optimal potential by extraction of the envelope of 
curves from {log (E% ^ (E v , 5, Ef n ))} E m in , calculated 
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log(E /GeV) 

FIG. 5: Power law limit representation of IceCube discovery 
potential. The thin curves (from left to right) corresponds to 
increasing reconstructed muon energy thresholds, from 100 
GeV to 100 TeV, in step of 0.5 in log(E), for a declination 
S = 15° . The thick curve is the envelope, representing the best 
possible reach following the selection method of reconstructed 
muon energy threshold. 



for the sources at different declinations 6. The mean- 
ing of each curve is the following: A power law neutrino 
spectrum with normalization tangent to the thick curve 
is at the discovery limit. 

Technically, the extraction of the envelope is per- 
formed in two steps instead of the single one as il- 
lustrated in Fig [5) once the best flux normaliza- 
tion $g un (r^ , S) determined (the lower envelope of 
{\og$o(T„, S, £™ in )}£mi„), the smooth thick curve from 
Fig [5] is obtained from the envelope extraction from the 
series of curves {log$Jf ln (r i/ ,(5)}r„. The envelope condi- 
tion reads 

H lno' <t> min 

log^^MF.,,5), (A2) 

dr„ 

thus 

log(E 2 u ^(E v ,S)) = log$™»(T„,5) + (2 - T v )logE„. 

(A3) 

The results are shown Fig. [4] for various declinations 
and represent the IceCube discovery potential to power 
law flux spectra. 
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